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This  is  the  first  of  a  series  of  papers  addressing  the  solution  of  the  inverse  scattering 
problem  for  the  Helmholtz  equation  in  two  dimensions.  Here,  we  derive  a  system 
of  differential  equations  for  the  scattering  matrices  which 

1.  Directly  govern  the  whole  behavior  of  the  scattering  problem, 

2.  Can  be  easily  implemented  numerically  with  any  prescribed  accuracy. 

In  the  subsequent  papers,  we  will  use  this  apparatus  to  design  stable  inversion  al¬ 
gorithms  for  the  acoustical  inverse  scattering  problem.  Specifically,  in  the  second 
paper,  we  will  present  a  scheme  based  on  the  trace  formula  which  is  a  direct  exten¬ 
sion  to  the  one  employed  in  [16].  The  algorithm  is  quite  satisfactory  analytically, 
but  requires  excessive  amounts  of  CPU  time.  Finally,  in  the  third  paper,  we  will 
present  a  radically  accelerated  version  of  the  algorithm. 
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1  Introduction 


During  the  last  several  decades,  the  inverse  scattering  problems  for  the  Helmholtz 
equation  have  enjoyed  a  remarkable  degree  of  popularity,  both  in  pure  and  applied 
contexts  (see,  e.g.,  [2],  [3]).  A  number  of  algorithms  have  been  proposed  for 
the  numerical  treatment  of  these  problems,  in  such  environments  as  medical 
diagnostics,  non- destructive  testing,  anti-submarine  warfare,  oil  exploration,  etc. 
In  the  design  of  such  a  scheme,  three  major  problems  have  to  be  overcome. 

1.  The  problem  is  highly  non-linear,  even  in  its  purely  mathematical  form.  In 
the  one- dimensional  case,  the  problem  can  be  reduced  to  a  linear  one,  but  the 
procedure  is  not  stable  numerically. 

2.  Once  a  mathematically  valid  inversion  scheme  is  constructed,  it  may  or  may 
not  be  stable  numerically.  In  fact,  no  numerically  robust  schemes  seem  to  exist 
at  this  time,  except  in  one  dimension. 

3.  The  cost  of  applying  the  scheme  on  the  computer  tends  to  be  extremely  high, 
except  in  the  one-dimensional  case. 

The  existing  attempts  to  solve  inverse  scattering  problems  for  the  Helmholtz 
equation  can  be  roughly  subdivided  into  four  groups. 

1.  Linearized  inversion  schemes,  attempting  to  approximate  the  inverse  scattering 
problem  by  the  problem  of  inverting  an  appropriately  chosen  linear  operator  (see, 
for  example,  [3]). 

2.  Methods  based  on  non-linear  optimization  techniques,  attempting  to  recover 
the  parameters  of  the  problem  iteratively,  by  solving  a  sequence  of  forward  scat¬ 
tering  problems  (see,  for  example,  [4],  [5],  [6]). 

3.  Gel’fand-Levitan  and  related  techniques,  converting  the  Helmholtz  equation 
into  the  Schrodinger  equation,  the  inverse  problem  for  the  latter  being  reducible 
to  the  solution  of  a  linear  Volterra  integral  equation  (see,  for  example,  [2],  [7]). 

4.  Techniques  based  on  the  trace  formulae,  connecting  the  high-frequency  be¬ 
havior  of  the  solutions  of  the  Helmholtz  equation  with  the  local  values  of  the 
parameters  to  be  recovered  (see,  for  example,  [10],  [11],  [12]). 

From  the  mathematical  viewpoint,  the  one- dimensional  problem  was  satisfac¬ 
torily  solved  in  the  early  fifties  (see  [7],  [8]).  However,  procedures  of  the  type 
described  in  [7]  and  [8]  do  not  lead  to  stable  numerical  algorithms.  Existing 
stable  and  efficient  schemes  in  one  dimension  are  based  on  the  trace  formulae. 
They  consist  of  constructing  a  Riccati  equation  for  some  function  of  the  solution 
of  the  Helmholtz  equation  (such  as  the  impedance  or  the  scattering  coefficient), 
and  combining  it  with  a  trace  formula.  The  result  is  a  system  of  ordinary  differ- 
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ential  equations  that  can  be  solved  numerically,  and  a  proper  choice  of  the  trace 
formula  ensures  stability  and  rapid  convergence  of  the  process  (see,  for  example, 
[10],  [11],  and  [16]). 

Attempts  to  generalize  this  approach  to  higher  dimensions  have  not  lead  to 
effective  numerical  procedures;  however,  a  collection  of  powerful  mathematical 
apparatus  has  been  developed  (see  [9]  and  [10]  for  quantum  scattering  in  one 
dimension;  [13],  [14]  and  [17]  for  electrical  impedance  scattering;  and  [15]  for 
Schrodinger  scattering  in  three  dimensions).  One  of  the  more  promising  tools 
developed  to-date  is  the  concept  of  the  scattering  matrix  (or  of  the  Dirichlet-to- 
Neumann  map)  and  the  differential  equations  governing  it  (see  [13],  [14],  [17]). 

The  present  work  is  the  beginning  of  a  series  of  papers  addressing  the  solution 
of  the  inverse  scattering  problem  for  the  Helmholtz  equation  in  two  dimensions. 
In  this  paper,  our  goal  is  to  derive  a  system  of  differential  equations  for  the 
scattering  matrices  which 

1.  Directly  governs  the  whole  behavior  of  the  scattering  problem, 

2.  Can  be  easily  implemented  numerically  with  any  prescribed  precision. 

In  the  subsequent  papers,  we  will  use  this  apparatus  to  design  the  inversion 
algorithms.  Specifically,  in  the  second  paper,  we  will  present  an  inverse  scattering 
scheme  based  on  the  trace  formula  which  is  a  direct  extension  to  the  one  employed 
in  [16].  The  algorithm  is  quite  satisfactory  analytically,  but  requires  excessive 
amounts  of  CPU  time.  Finally,  in  the  third  paper  to  be  published,  we  will  present 
a  much  accelerated  version  of  the  algorithm. 

The  principal  results  of  the  paper  are  the  Riccati  equations  for  the  scattering 
matrices  in  cylindrical  coordinates.  The  paper  is  organized  as  follows.  In  Sec¬ 
tion  2,  we  summarize  the  relevant  properties  of  the  Helmholtz  equation  in  two 
dimensions,  and  introduce  the  scattering  matrices.  In  Section  3,  we  derive  the 
differential  equations  which  the  scattering  matrices  satisfy. 


2  Mathematical  Preliminaries 


In  this  section,  we  will  discuss  the  Helmholtz  equation  and  its  associated  scatter¬ 
ing  problems.  First  in  Section  2.1,  we  introduce  common  and  special  usages  of 
notation  in  this  paper.  In  Section  2.2,  properties  of  the  Bessel  functions  are  pre¬ 
sented.  Section  2.3  is  then  devoted  to  the  scattering  problems  for  the  Helmholtz 
equation.  Finally  in  Sections  2.5  and  2.6,  we  define  scattering  matrices  corre¬ 
sponding  to  three  special  scattering  problems. 
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2.1  Notation 


We  will  denote  by  C'*"  the  set  of  all  complex  numbers  with  nonnegative  imaginary 
part,  and  by  5^  C  C  the  unit  circle  in  the  complex  plane  defined  by  the  formula 

I  2:  1=  1;  (1) 


we  will  assume  that  5"^  is  parameterized  by  its  arc  length.  We  will  define  i  as 
the  linear  space  of  all  two-sided  sequences  of  complex  numbers 

^  =  {  =  0,±1,±2,- •  •  }.  (2) 

In  agreement  with  standard  practice,  we  will  denote  by  the  subspace  of  £ 
consisting  of  all  sequences  ^  such  that 

00 

Y1  I  l^<  (3) 

m=“Oo 

and  by  £°°  the  subspace  of  £  consisting  of  all  sequences  ^  such  that 

sup  I  (rn  |<  00.  (4) 


Let  F  :  L'^{S^)  £'^  denote  the  Fourier  transform  converting  a  square 

integrable  function  on  the  circle  into  its  Fourier  series,  so  that  the  expression 

+  00 

m=  E  f".''""  (5) 

m=  — 00 

can  be  written  in  the  matrix  form 

/  =  (6) 

When  it  is  necessary  to  explicitly  show  the  dependence  of  on  the  variable  6, 
the  subscript  0  will  be  suffixed  to  the  linear  mapping  F~^-,  therefore  the  expression 
(6)  can  be  rewritten  as 

m = (7) 

In  agreement  with  standard  practice,  we  denote  by  Jm  the  Bessel  function 
of  the  first  kind  of  order  m,  and  by  Hm  the  Hankel  function  of  the  first  kind  of 
order  m.  We  will  denote  by  J^,  Hz  the  infinite  diagonal  matrices 

Jz  =  diag{...,J-i{z),jQ{z),Ji{z),...},  (8) 

Hz  =  diag{...,H.^{z),Ho{z),H^{z),...}.  (9) 
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Given  i?  >  0,  we  will  denote  by  D{R)  the  disk 

D{R)  =  {{r,e)  I  t<R},  (10) 

by  E{R)  the  exterior  of  D{R) 

E{R)  =  {  (r,^)  I  r>R  },  (11) 

and  by  A{R,  h)  for  h  >  0  the  annulus 

A[R,  /i)  =  {(r,  0)  I  +  (12) 

Further,  for  an  arbitrary  a:  e  we  will  denote  by  D{x,  R)  the  disk  of  radius  R 
centered  at  x.  For  a  continuous  function  q  €  C{D{R)),  and  a  real  number  r  <  R, 
we  define  the  mapping  :  C{S^)  C(5^)  by  the  formula 

iQr-fm^q{r,0)-f{9).  (13) 

Given  a  function  q:  R'^  R^,  and  a  set  C  we  will  define  the  function 

the  restriction  of  q  on  A,  via  the  formula 

Qa{x)  =  q{x)  ■  xiA),  (14) 

with  x(>l)  the  characteristic  function  of  A. 

For  an  arbitrary  z  ^  C,  we  will  denote  by  the  linear  space  of  all  two-sided 
complex  sequences  {l^m}  such  that  for  some  c  >  0, 


for  all  integer  m.  We  will  denote  by  the  linear  space  of  all  two-sided  complex 
sequences  {am]  such  that  for  some  c  >  0, 


for  all  integer  m. 

2.2  Several  Classical  Lemmas 

In  this  subsection,  we  summarize  several  classical  results  describing  the  behavior 
of  certain  special  solutions  of  the  Helmholtz  equation.  The  following  lemma  can 
be  found,  for  example,  in  [21]. 
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Lemma  2.1  Let  m  be  an  integer.  Suppose  that  Jm.  is  the  Bessel  function  of 
order  m,  and  Ym  is  Neumann  function  of  order  m.  Suppose  further  that  Hm  is 
the  first  kind  Hankel  function  of  order  m  defined  by  the  formula 

=  +  Ym.  (17) 


Then 

(I)  Jm  =  (-1)"^^-^,  and  Hm  =  {-l)”^H.m. 

(II)  For  any  z  €  C'"*', 

with  7  =  0.5772 . . .  the  Euler  constant. 

(III)  For  any  z  G  , 

_  /  2?t 

lim  Jmiz)  •  V2nm  •  (  — 

m—^oo  \  6  *  . 

and 

Hm  Hm{z)-^-(^ 

m-*oo  ^  \  2  \  2n 

In  other  words,  for  a  fixed  z  G  as  m  increases,  Jm{z)  decays  as 

V^Trm  \ 

and  Hm{z)  grows  as 


2m 


=  1, 


Hm{z)^-iJ—-(—)  . 

V  Trm  \e  •  zj 

(IV)  For  a  fixed  m  and  large  z  G  C'^ ^ 


^rn{z)  =  \l^{^^v[i-{^-m^-l)]+0{z-^)] 

More  specifically,  for  z  G  C'^  and  \z\  >  2, 


(18) 

(19) 


(20) 

(21) 

(22) 

(23) 

(24) 

(25) 

(26) 
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with 


1 


(27) 


4^)\  < 


8  •  \z 


The  following  lemma  provides  an  upper  bound  for  the  Hankel  function  Hq. 


Lemma  2.2  Suppose  that  z  G  is  a  complex  number  in  the  upper  half  of  the 
complex  plane.  Then 

I  ^o(^)  |<  2  (28) 


for  all  \  z  \>  2,  and 


for  all  \  z  \<  2. 


Ho{z)  \<  2+  1  ln(|  z/2  |)  | 


(29) 


Proof.  For  \z\  >  2,  formula  (28)  follows  directly  from  (26).  For  1^1  <  2,  formula 
(29)  follows  directly  from  the  combination  of  (17),  (18)  and  (19).  □ 

For  a  domain  0,  C  a,  point  x  £  R?,  and  a  positive  number  k,  we  will  denote 
by  |f2|  the  area  of  fl,  and  by  P{x,  0)  the  real  number  defined  by  the  formula 

P{x,  ^)=  I  I  ln(^lk  -  II)  I  (30) 

JQnD(x,l/k) 

The  following  technical  lemmas  provide  estimates  for  P(x,Q). 


Lemma  2.3  For  arbitrary  domain  fl  C  R^ ,  point  x  G  Rf ,  and  positive  number 
k  G  R\ 

P{x.,D,)  <  P{x,D{x,p))  (31) 


provided  that 


P  = 


(32) 


Proof.  It  obviously  follow  from  (32)  that 

1.  The  domain  0.  and  the  disk  D{x,p)  have  the  same  area, 

2.  D{x,p)  is  contained  in  D{x,llk) 

Assuming  that  0  does  not  coincide  with  D{x,p),  that  is,  assuming  that  the  two 
domains 


(33) 

(34) 
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A  =  n\D(x,p}, 

B  =  p) 


have  non-zero  areas  |yl|  =  \B\,  we  intend  to  show  that 


i\n\ 


(35) 


Without  loss  of  generality,  we  may  assume  that  Q,  is  contained  in  D(x,  1/k).  It 
follows  from  (32)  that  for  any  (fi  €  A  C  D{x,  1/k)  and  ^2  €  S  C  D{x,  1/k) 


we  have 
Therefore, 


k\\x  -  I2II  <k-  p<  -  6||  <  1, 

Ik -611  <  Ik -611- 


I  M^lk-611)  l<l  M^lk-611)  I 

which  establishes  (31).  □ 

Lemma  2.4  For  arbitrary  domain  Q,  C  and  point  x  G  R^, 


P(x,f))<i^|l  + 


In  I  I 


TT 


(36) 

(37) 

(38) 


(39) 


Proof.  We  first  consider  the  case  when  the  area  of  O  is  no  less  than  that  of  the 
disk  D{x,l/k),  that  is, 

l^^l  >  (40) 

Under  this  assumption,  it  immediately  follows  from  (30)  that 


P{x,  fl) 


< 


L 


DiT,l/k) 

27r  fl/k 


ln(fc|k  -61/2)  I 


ni/K 

\n{kr)rdrd9 

_ 


JL  <  M 

2P  -  2  ’ 


(41) 


from  which  (39)  follows.  Now,  we  consider  the  case  when  the  area  of  fi  is  less 
than  that  of  the  disk  D{x,  1/k),  that  is, 


Let 


(42) 

(43) 
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According  to  Lemma  2.3, 


□ 


P{x,fl)  < 


P{x,  D{x,  p))  =  f  I  \n{k\\x  -  ^11)  I 

JD(x,p) 

nP 

\n(kr)rdrd9 

. 


M 

2 


1  +  2- 


(44) 


Definition  2.5  In  agreement  with  standard  usage,  we  will  refer  to  the  function 
f  :  C  in  the  Helmholtz  equation 

A<f>  +  e(f>  =  f  (45) 

as  the  source,  and  to  the  solution  cf)  :  Ri^  ^  C  as  the  field  generated  by  the  source. 


The  following  two  lemmas  can  be  found,  for  example,  in  [1]. 

Lemma  2.6  Suppose  that  k  G  C'^  is  a  complex  number.  Then  the  Greenes  func¬ 
tion  for  the  homogeneous  Helmholtz  equation 

+  ArV  =  0  (46) 

subject  to  the  Sommerfeld  radiation  condition 

=0  (47) 

is  given  by  the  formula 

Gi(x,0  =  -j/7o(*||i-e||),  (48) 

where  x  =  {r,6)  and  ^  =  {r',9')  are  arbitrary  distinct  points  in  . 


Lemma  2.7  (Graf’s  Addition  Formula)  Suppose  that  the  positive  numbers  u,v, 
w,  a  are  such  that  u  >  v,  and  that 


—  2uv  cos(q:). 

Then 

CO 

Ho{w)  =  Hm{u)Jm{v)e'”’''’‘ 

7n=  — CO 


(49) 

(50) 
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The  following  lemma  is  an  immediate  consequence  of  Lemma  2.7. 


Lemma  2.8  Suppose  that  x,  ^  are  two  points  in  such  that  ||x||  >  ||^||.  Sup¬ 
pose  further  that  {r,9)  and  {r',9')  are  the  polar  coordinates  of  the  vectors  x, 
respectively.  Then 


OO 

^^o(*:||x-{||)=  E  H„(kr)Mkr')e‘’'^’-‘'K  (51) 

m=— OO 


The  following  lemma  is  a  direct  result  of  Lemmas  2.6  and  2,2. 


Lemma  2.9  Suppose  that  k  G  is  non-zero,  and  that  ft  G  B?  is  a  domain  of 
area  p.  Then 


f  I  Gi(x,{)  I  di 
Jq 


< 


2  “ 


(52) 


for  all  X  ^  B? . 


Proof.  Introducing  the  notation 

I(x)=  j  \G,(x,(,)\di  (53) 

and  using  (48),  we  have 

I{x)  =  Uf  \  H„{z)  \  d(  +  ^  \  H„{z)  \  dA  ,  (54) 

with  ^  =  is;||x  -  ^11/2.  Combining  (28)  and  (29)  with  (54),  we  immediately  obtain 
the  estimate 

l{x)  <  1  (2^  +  /  I  lii(U  I)  i  dA  (55) 

4  \  JQn{\z\<i}  J 

for  all  X  ^  R? .  Now  the  lemma  follows  immediately  from  (55)  and  Lemma  2.4. 
□ 


Remark  2.10  Denoting  by  Gk,Q  the  linear  operator  (7(11) 
the  the  formula 

{Gk,Q-^){x)=  [  Gk{x,()'tp(t)d^, 
jQ 

for  all  Ip  G  C{D,),  we  can  rewrite  (52)  in  the  form 


lia 


A:,0l|oo  — 


C{0,)  defined  by 

(56) 

(57) 
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The  following  lemma  provides  a  sharper  upper  bound  for  ||GA:,n||oo  when  is  the 
annulus  A{R,h)  (see  (12)).  Its  proof  is  nearly  identical  to  that  of  Lemma  2.9 
and  is  omitted. 

Lemma  2.11  Suppose  that  i?  >  0,  h  >  0  are  two  real  numbers,  and  that  k  G  C'*' 
is  non-zero.  Then  for  all  x  ^  B? 

I (58) 

where  fi  =  h  •  ir  •  {2R  +  h)  is  the  area  of  the  annulus  A(i?,  /i). 

Remark  2.12  The  preceding  two  lemmas  show  that  there  are  two  types  of  esti¬ 
mates  (see  formulae  (52)  and  (58))  on  ||G/:,q||oo;  depending  on  the  shape  of  the 
domain.  When  the  domain  is  a  disk,  the  estimate  is  of  the  form 

||G^,fi||co  <  c-  |ln(//)|  •  (59) 

whereas  when  the  domain  a  annulus,  the  estimate  is  of  the  form 

||G^/:,q||co  <  c  -  p,  (60) 

Definition  2.13  A  function  (j)  is  said  to  be  a  radiation  field  in  a  bounded  domain 
Cl  if  and  only  if  (f)  is  a  solution  of  the  homogeneous  Helmholtz  equation  (f6)  in 
Cl;  a  function  cf)  is  said  to  be  a  radiation  field  outside  a  bounded  domain  Cl  if  and 
only  if  ({)  is  a  solution  of  the  homogeneous  Helmholtz  equation  (f6)  subject  to  the 
Sommerfeld  radiation  condition  (47), 

Lemma  2.14  Suppose  that  k  is  an  arbitrary  complex  number,  R  is  a  positive 
real  number,  and  (j)  :  D{R)  C  is  a  radiation  field  in  D{R).  Then  there  exists 
a  sequence  of  complex  numbers  {  a;  },  j  =  0,  ±1,  ±2, . .  such  that 

oo 

(f>{r,e)=  c^mMkr)e'^^  =  (61) 

m=  — OO 

for  all  r  <  R  and  0  <  0  <  27r. 

The  following  two  lemmas  are  widely  known.  They  can  be  found,  for  example, 
in  [1], 

Lemma  2.15  Suppose  that  k  is  an  arbitrary  complex  number,  R  is  a  positive 
real  number,  and  (f)  :  R'^  \  D{R)  C  is  a  radiation  field  outside  D(R).  Then 
there  exists  a  sequence  of  complex  numbers  {  aj  },  j  =  0,  ±1,  ±2, . . .,  such  that 

OO 

^(r,0)=  ^  /3mHm{kr)e'”^^  =  Fi^Hkr/d,  (62) 

m=— OO 

for  all  r  <  R  and  0  <  0  <  27r. 
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Remark  2.16  It  immediately  follows  from  Lemma  2.1  that  the  series  (61)  con¬ 
verges  inside  D{R)  if  and  only  if  a  E  YkR,  and  the  series  (62)  converges  outside 
D{R)  if  and  only  if  ^  E  XkR. 

2.3  Scattering  Problems 

The  subject  of  this  paper  is  the  Helmholtz  equation 

+  g)<^  =  0.  (63) 

In  (63),  we  assume  that  A;  is  a  complex  number  such  that  Im{k)  >  0,  and 
q  :  R^  R^  is  a.  smooth  function;  we  will  be  referring  to  the  function  q  as  the 
scatterer.  We  further  assume  that  the  support  of  g  is  a  bounded  domain  Q,  and 
that  q{x)  >  —  1  for  all  x  G  R^.  We  will  be  considering  solutions  H  i— >•  C  of 
Equation  (63)  of  the  form 

=  <l>o  -\-  '4’ 1  (64) 

with  4>q  a  radiation  field  in  fl,  and  4)  :  R?  ^  C  a  radiation  field  outside  fl.  We 
will  be  referring  to  as  the  total  field,  to  4>o  as  the  incoming  field,  and  to  ^  as 
the  scattered  field.  Furthermore,  we  will  be  referring  to  the  determination  of  the 
scattered  field  from  a  given  incoming  field  as  the  (forward)  scattering  problem. 
It  is  easy  to  verify  that  tp{x)  satisfies  the  equation 

A4>  +  k'^4’  =  —k^Qi4>o  +  4’)  (65) 


for  all  X  E  R?. 

Remark  2.17  It  is  well-known  (see,  for  example,  [19],  [18],  [20])  that  the  for¬ 
ward  scattering  problem  is  well-posed.  More  specifically,  the  problem  can  be  re¬ 
formulated  as  the  so-called  Lippmann-Schwinger  equation 

4^{x)  +  e  f  Gk{x,0qm{0d^  =  -k^  I  (66) 

JQ  Jn 

for  the  scattered  fields  or  as  the  Lippmann-Schwinger  equation 

4>{x)  +  k^  I  Gk{x,  0q{04>{0d^  =  <l>o{x).  (67) 

J 

for  the  total  field,  for  all  x  E  R^;  either  of  the  second  kind  integral  operators  (66), 
(61)  is  invertible,  and  the  maximum  norm  of  the  inverse  operator  is  bounded. 
That  is  to  say,  for  k  E  and  q  E  C{R?)  having  the  compact  support  ft,  an 
incoming  field  ([q  determines  uniquely  a  scattered  field  xp,  and  xp  depends  contin¬ 
uously  on  <po  in  maximum  norm. 
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Defining  the  linear  operator  Gl  :  C{Q.)  *->■  C{B?)  by  the  formula 


(Gl  ■  i,)(x)  =  e  [  Gi(x,e)?(OV’K)<if, 

(68) 

we  can  rewrite  (66)  in  the  form 

V’  +  Gli;  =  g, 

(69) 

with  g  €  C{B?)  defined  by  the  formula 

g  =  G\fQ. 

(70) 

When 

l|Gll|oo<l, 

the  equation  (69)  can  be  solved  via  the  fixed-point  iteration 

(71) 

V’o  =  0, 

i’m+l  =  g  -  Gld>m- 

(72) 

(73) 

It  follows  from  Lemma  2.9  that  the  condition  (71)  is  met  whenever  the  area  of 
the  domain  fl  is  sufficiently  small.  The  above  discussion  is  formalized  in  the 
following  lemma  whose  simple  proof  we  omit. 

Lemma  2.18  Suppose  that  k  €  C"*"  is  a  complex  number,  q 
is  a  domain  of  area  p,.  Suppose  further  that 

€  C{Rf),  and  Cl  C 

f -f  llnfl  yt  1  •,/^)| 

=  — ' 4 

(74) 

Finally,  suppose  that  p  is  sufficiently  small  so  that 

8<\. 

(75) 

Then  the  fixed-point  iteration  (72),  (73)  converges  to  the  solution  of  Lippmann- 
Schwinger  equation  (69).  Moreover, 

\  '^{x)  -  V’m(a;)  |<  6”"+^ 

(76) 

for  all  integer  m  >  0,  and  x  ^  B? . 


When  fl  is  an  annulus  A(R,h)  with  h  small,  Lemma  2.18  assumes  the  form 
provided  by  the  following  lemma. 
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Lemma  2.19  Suppose  that  R  >  0  and  h  >  0  are  two  real  numbers.  Suppose 
further  that  k  €  is  a  complex  number,  q  G  C{R?),  and  Q  =  A(R,h)  is  an 
annulus  of  width  h.  Finally,  suppose  that  h  is  less  that  1,  and  that 


|tn|,U2B+l)’ 

Then  the  fixed-point  iteration  (72),  (73)  converges  to  the  solution  of  Lippmann- 
Schwinger  equation  (69).  Moreover, 

I  'Ipix)  -  tpmix)  |<  (tt-  I  k  p  •||9||oo(2i?  +  1)  •  (78) 

for  all  integer  m  >0,  and  x  ^  R? . 


2.4  Scattering  in  Circular  Geometry 

In  the  remainder  of  this  paper,  we  will  be  interested  in  two  special  cases  of  the 
scattering  problem:  in  the  first  case,  the  compact  support  of  the  scatterer  q  is 
a  disk  D{p)  with  some  p  >  0;  in  the  second  case,  the  compact  support  of  the 
scatterer  is  a  domain  ft  lying  outside  the  disk  D(p),  that  is,  Cl  C  R?  \  D{p),  see 
Figure  1. 

In  this  subsection,  we  construct  simple  analytical  expressions  for  the  incoming 
and  the  scattered  fields  in  each  of  these  two  cases.  Clearly,  in  the  first  case,  the 
incoming  field  (f>o  is  a  radiation  field  inside  D{p),  and  is  therefore  generated  by 
sources  located  outside  D{p)-,  the  scattered  field  ^  is  a  radiation  field  outside 
D{p),  and  is  therefore  generated  by  sources  inside  D{p).  We  will  be  referring  to 
this  scattering  problem  as  the  interior  scattering  problem  (see  Figure  l:(a)).  The 
following  lemma  is  a  direct  consequence  of  Lemmas  2.14,  2.15. 

Lemma  2.20  Suppose  that  p  is  a  positive  real  number,  and  that  the  scatterer  q 
has  the  compact  support  D{p).  Then 

(I)  If  the  function  <j>o  :  D(p)  C  is  an  incoming  field  to  the  scatterer  q,  then 
there  exists  a  sequence  a  G  Ykp,  such  that  for  all  {r,9)  G  D{p), 

OO 

amJm{kr)e''^^  =  Ff'^JkrOf,  (79) 

m=  — OO 

(II)  If  the  function  xj)  •.  Rf  C  is  a  scattered  field  from  the  scatterer  q,  then 
there  exists  a  sequence  fl  G  Xkp,  such  that  for  all  {r,$)  G  \  T){p), 

OO 

xP{r,e)=  (3mHm{kr)e'^^  =  Fg-^HkrlS.  (80) 

m=— OO 
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Figure  1:  (a)  Interior  scattering  where  scatterer  occupies  the  entire  disk  D{p). 
(b)  Exterior  scattering  where  the  scatterer  is  located  outside  the  disk  D{p). 


In  the  second  case  of  interest  to  us,  we  consider  incoming  field  generated  inside 
the  disk  D{p),  and  therefore  the  incoming  field  is  a  radiation  field  outside  the  disk. 
Since  the  scattered  field  ^  is  a  radiation  field  outside  17,  it  is  a  radiation  field  inside 
D{p).  We  will  be  referring  to  this  scattering  problem  as  the  exterior  scattering 
problem  (see  Figure  l:(b)).  The  following  lemma  is  a  direct  consequence  of 
Lemmas  2.14,  2.15. 

Lemma  2.21  Suppose  that  p  is  a  positive  real  number,  and  that  the  scatterer  q 
has  a  compact  support  outside  D{p).  Then 

(I)  If  the  function  <j>o  :  Bf  \  D{p)  i— >  C  is  an  incoming  field  to  the  scatterer  q, 
then  there  exists  a  sequence  /3  G  Xkp,  such  that  for  all  {r,0)  G  \  D{p), 

OO 

Mr',0)=  E  /3mHrnikr)e'’^^  =  Fi^Hkr/3-,.  (81) 

m=-oo 


(II)  If  the  function  rp  •.  C  is  a  scattered  field  of  the  exterior  scattering 

problem,  then  there  exists  a  sequence  a  G  Ykp,  such  that  for  all  {r,0)  G  D{p), 

OO 

Y.  OimJm{kr)e"^^  =  Ff^JkrOc-  (82) 

m=  — OO 

2.5  Scattering  Matrices 

One  of  principal  analytical  tools  used  in  this  paper  are  the  scattering  matrices; 
a  scattering  matrix  is  the  linear  mapping  converting  the  incoming  field  into  the 
scattered  field,  for  a  specified  scatterer.  For  technical  reasons,  it  is  convenient  to 
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have  slightly  different  definitions  of  the  scattering  matrices  for  the  interior  and 
exterior  scattering  problems;  thus,  we  define  the  interior  and  exterior  scattering 
matrices  separately. 

Definition  2.22  Under  the  conditions  of  Lemma  2.20,  any  a  G  YkR^  defines  an 
incoming  field  (f>o  via  (79);  the  resulting  scattered  field  i/)  can  be  represented  in  the 
form  (80),  with  €  XkR^.  Thus,  there  exists  a  linear  mapping  Si  :  YkR^  XkFu, 
such  that 

=  (83) 

for  all  a  G  YkR^)  and  we  will  refer  to  SI  as  the  interior  scattering  matrix  corre¬ 
sponding  to  the  scatterer  q  and  frequency  k. 

Remark  2.23  Obviously^  for  a  fixed  k  £  C'^ j  the  scattering  matrix  contains  all 
the  information  that  can  be  acquired  by  any  scattering  experiments  performed 
outside  the  scatterer;  we  will  refer  to  the  collection  of  the  scattering  matrices 

{  I  kiR')  (84) 


as  the  complete  scattering  data. 

Remark  2.24  In  an  actual  scattering  experiment^  measurements  are  obtained 
at  a  finite  collection  of  frequencies  kj.  Therefore,  in  a  more  realistic  formula¬ 
tion  of  the  inverse  scattering  problem,  the  scatterer  q  is  to  be  determined  from 
measurements  of  the  scattering  matrices  at  finite  number  of  frequencies: 

{SJJ  j  =  1.2....,Ar}.  (85) 

In  the  remainder  of  this  paper,  we  will  assume  that  our  scatterer  q  has  compact 
support  D{Ro)  for  some  positive  number  Rq,  so  that 

9(r,0)  =  O  (86) 

for  all  r  >  i?o-  Furthermore,  for  a  positive  number  R,  we  will  denote  by 
the  interior  scattering  matrix  corresponding  to  the  scatterer  qD{R)  (see  (14))  and 
frequency  k\  in  other  words, 

s-R,,  =  sr""-  m 

We  will  refer  to  the  function  qr)(R)  as  the  truncated  scatterer.  Obviously,  at 

i?  =  0,  the  truncated  scatterer  is  zero;  it  therefore  generates  no  scattered  field; 
in  other  words, 

SZ.k  =  0.  (88) 

Likewise,  since  qD(R)  =  Q  for  any  R  >  Ro,  we  have 

^R,k  —  (S9) 


15 


for  all  R  >  Rq. 


Closely  related  to  operators  Sp^  f.  are  the  so-called  exterior  scattering  matrices 
associated  with  the  scatterer  qE{R)  (see  (11), (14)). 

Definition  2.25  Under  the  conditions  of  Lemma  2.21,  any  jd  G  XkR  defines  an 
incoming  field  4>o  via  (81);  the  resulting  scattered  field  if  can  be  represented  in 
the  form  (82),  with  some  a  €  YkR-  Thus,  there  exists  a  linear  mapping  : 
YkR  i-->  XkR,  such  that 

(90) 

for  all  (3  G  and  we  will  refer  to  as  the  exterior  scattering  matrix  corre¬ 

sponding  to  the  scatterer  qE(R)  o,7id  frequency  k. 

Remark  2.26  We  will  refer  to  qE(R)  cls  the  hollowed  scatterer.  Since  the  scat¬ 
terer  q  has  compact  support  in  D{Ro)^  the  hollowed  scatterer  qE{R)  is  zero  outside 
D{Ro).  Therefore^ 

Stk  =  0  (91) 

for  all  R  >  Rq. 


2.6  Scattering  from  an  Annulus 

Given  a  pair  of  real  numbers  R  >  0  and  h  >  0,  we  will  refer  to  a  scatterer 
q  :  Rf  t—y  R^  a.s  an  annular  scatterer  of  inner  radius  R  and  width  h  if  q{r,6)  —  0 
for  all  r  <  i?  or  r  >  i?  +  h.  Obviously,  for  any  scatterer  q,  the  scatterer  qA{R,h) 
(see  (12), (14))  is  an  annular  scatterer. 

Conceptually,  the  incoming  field  to  the  annular  scatterer  qA(R,h)  is  generated 
by  sources  both  inside  the  disk  D{R)  and  outside  the  disk  D{R  -|-  h)  since  it  is  a 
radiation  field  inside  the  annulus  A{R,  h).  By  the  same  token,  the  scattered  field 
from  the  annular  scatterer  qA(R,h)  is  a  radiation  field  outside  the  annulus  A(R,  h), 
and  therefore  has  the  form  (61)  inside  the  disk  D{R),  and  the  form  (62)  outside 
the  disk  D{R  +  h).  The  following  two  obvious  lemmas  formalize  these  facts. 

Lemma  2.27  Suppose  that  R,  h  are  two  positive  numbers,  and  that  A{R,h)  is 
the  annulus  defined  via  (12).  Then  (f>o  :  A{R,h)  C  is  an  incoming  field  to 
A{R,  h)  if  and  only  if  there  exist  two  functions  :  D{R  +  h)  C,  : 
R^  \  D{R)  I— »  C  such  that 

=  +  (92) 

where  is  a  radiation  field  in  D{R  +  h);  in  other  words,  there  uniquely  exists 
O!  €  Y^Rph'j,  such  that  for  all  {r,9)  G  D{R  +  h), 

,^(-)(r,0)  =  F,-Vfc.o;  (93) 
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and  where  is  a  radiation  field  outside  D{R);  in  other  words,  there  uniquely 
exists  (3  G  XkR,  such  that  for  all  {r,6)  €  \  D{R), 

=  Ff^Hkrfi.  (94) 

Furthermore^  the  decomposition  (92)  is  unique. 

Lemma  2.28  Under  the  conditions  of  Lemma  2,27  Suppose  that  'ip  :  ^  C  is 

the  scattered  field  generated  by  the  annular  scatterer  A(R^h),  Then  there  exist 
two  sequences  6l  G  YkR  and  G  such  that  inside  the  disk  D{K)j 

xlf{r,e)  =  Ff^JkrC^,  (95) 

and  outside  the  disk  D{R  +  h), 

^{r,e)  =  F,-'H„ff.  (96) 

Remark  2.29  Abusing  the  notation  somewhat,  we  will  be  referring  to  the  coef¬ 
ficients  a,fi  in  (93),  (94)  as  the  incoming  potential,  the  coefficients  a,$  in  (95), 
(96)  as  the  outgoing  potential. 

Obviously,  given  an  annular  scatterer  qA(R,h)i  the  combination  of  Lemmas  2.27 
and  2.28  defines  a  linear  mapping  Lfc(K+/i)  x  XkR  1-^  Xk(R+h)  x  YkR  converting 
the  incoming  potential  (a,  fi)  G  YkR  x  XkR  into  the  outgoing  potential  (fi,  a)  G 
XkR  X  YkR.  Thus,  we  are  led  to  the  following  definition. 

Definition  2.30  For  an  annular  scatterer  qA(R,h)!  scattering  matrix  S^  k 
the  mapping  Yk(R+h)  x  XkR  ^  Xk(R+h)  x  YkR  such  that  for  any  incoming  potential 
the  outgoing  potential  is  given  by  the  formula 


(97) 


Remark  2.31 

submatrices 


so  that 


The  scattering  matrix  can  obviously  be  partitioned  into  four 


skk  = 

Saa  SccfS 

(98) 

SjSaOf  +  Si3f3fi,) 

(99) 

a  = 

SoeO(^  “t”  Sapfi- 

(100) 

17 

2.7  Scattering  from  Two  Disjoint  Scatterers 

When  several  disjoint  scatterers  are  irradiated  by  an  incoming  field,  we  wish  to 
calculate  the  scattered  field  under  the  condition  that  we  know  how  to  solve  the 
scattering  problem  associated  with  each  scatterer  separately.  In  other  words, 
if  we  are  given  the  scattering  matrices  of  individual  scatterers,  we  would  like 
to  combine  them  and  obtain  the  scattering  matrix  associated  with  the  whole 
ensemble  of  the  scatterers.  In  this  subsection,  we  build  the  necessary  analytical 
machinery.  Although  the  discussion  is  confined  to  two  disjoint  scatterers,  all 
conclusions  can  easily  be  extended  to  the  case  of  multiple  scatterers. 

2.7.1  Lippmann-Schwinger  equation  for  a  Part  of  the  Scatterer 

Suppose  that  the  scatterer  q  consists  of  two  separate  scatterers  occupying  two 
disjoint  domains  fli  and  D,2: 


Oinflj  =  0,  (101) 

fli  U  ris  =  fl.  (102) 

Then  the  Lippmann-Schwinger  equation  (67)  can  be  rewritten  either  as 

+  k^  (  Gk{x,  i)q{i)mdi  =  M^)  I  Gk{x,  (103) 

Jill  Jil2 

or  as 

(l){x)  +  [  Gk{x,^)q{()<j){^)d^  =  (f>o{x)  -  k^  f  Gk{x,()q{^)<l){^)d^.  (104) 

JQi 

Introducing  the  functions  xj^i  :  B?  C  and  'ip2  '  ^  C  via  the  formulae 

xpi{x)  =  -k^  f  Gkix,0q{()4>{()d^,  (105) 

J 

x/>2{x)  =  -k^  [  Gk{x,C)qiO<l>U)d^,  (106) 

J  ii2 

we  observe  that  (103)  can  be  viewed  as  the  Lippmann-Schwinger  equation  (see 
(67))  on  the  domain  fli  for  the  total  field  <f>  corresponding  to  the  incoming  field 

<?^01  =  ^0  +  V’2-  (lOf) 

Similarly,  we  observe  that  (104)  can  be  viewed  as  the  Lippmann-Schwinger  equa¬ 
tion  on  the  domain  ^2  for  the  total  field  cf)  corresponding  to  the  incoming  field 

<?^02  =  ?!>0  +  V’l-  (108) 

The  following  obvious  lemma  formalizes  the  above  discussion. 
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Lemma  2.32  Suppose  that  the  scatterer  q  consists  of  two  separate  scatterers 
occupying  two  disjoint  domains  and  0,2-  Suppose  further  that  j>o  :  Q,i  U  Q2 
C  is  the  incoming  field  to  the  scatterer  q  and  that  <f>  :  Q,i  U  0.2  C  is  the 
corresponding  total  field.  Finally,  suppose  that  the  four  functions  ifj  :  Ef  C , 
<l)oj  Oj  I— C ,  for  j  =  1,2,  are  defined  by  the  formulae  (105)-(108).  Then 

(I)  For  j  =  1,2,  is  the  scattered  field  corresponding  to  the  incoming  field 
<f>Qj  to  the  scatterer  in  Oj. 

(II)  The  incoming  field  to  one  of  the  scatterers  is  the  superposition  of  the 
original  incoming  field  <f>Q  and  the  scattered  field  from  the  other  scatterer: 

<^01  =  ^0  +  V’2)  (109) 

<^02  =  +  (110) 

(III)  The  total  field  is  the  superposition  of  the  original  incoming  field  and  the 
scattered  fields  from  the  two  disjoint  scatterers: 

^  +  V’l  +  V’2-  (111) 

2.7.2  Merging  Two  Scatterers 

In  this  subsection,  we  examine  the  process  by  which  the  Lippmann- Schwinger 
equations  for  the  regions  Oi  and  O2  are  merged,  producing  the  Lippmann- 
Schwinger  equation  for  the  domain  Oi  U  172*  It  turns  out  that  when  the  area 
of  fli  is  small,  the  result  assumes  a  particularly  simple  form  (see  Lemmas  2.34 
and  2.37). 

Defining  the  linear  operators  Gj  :  C{Oj)  1— »■  C{Bf)  by  the  formula 

{Gr<i’)(x)  =  h'‘ !  (112) 

for  j  =  1,2,  and  using  the  formula  (111),  we  rewrite  the  Lippmann-Schwinger 
equations  (103),  (104)  in  the  form 

(7+G'i)-^i  +  G'i-V’2  =  (113) 

^2  •  +  (7  +  G2)  •  l/’2  =  —G2  •  <f>0-  (114) 

The  following  lemma  is  an  immediate  consequence  of  Remark  2.17. 

Lemma  2.33  Under  the  conditions  of  Lemma  2.32,  suppose  that  the  operators 
Pj  :  C{Ilj)  1-^  CR^),  Bj  :  C{Q,j)  1— >  C{R^)  are  defined  by  the  formulae 

P,  =  (/  +  G,)-',  (115) 

Bj  =  -(/+G,)-'-G,  =  -P,-G,,  (116) 

for  j  =  1,2,  Then  each  of  the  operators  Pi,  P2,  Pi,  P2  bounded  in  the 
maximum  norm. 
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Lemma  2.34  Under  the  conditions  of  Lemma  2.32,  suppose  that  the  estimate 
(60)  is  valid  on  Qi;  that  is,  there  exists  a  constant  c  such  that 

||G;t.nJ|oo  <  c-^,  (117) 

where  p  is  the  area  of  Lli.  Then,  for  all  sufficiently  small  p,  there  exist  unique 
functions  ^  C  and  ^  C  such  that  (113)  and  (II4)  are  satisfied. 

Furthermore, 

'ipi{x)  =  Bi  •  (I B2)(f>o -\- 0{p^),  (118) 

=  B2- {  I  +  Bi  ■  {I  +  B2)  }(l>o  +  0{p^),  (119) 

+  =  {  B2  +  iI  +  B2)-Br-{I  +  B2)}<l>o  +  0{p^).  (120) 

Proof.  Rewriting  (113)  and  (114)  in  the  form 

V’l  =  Bi{xp2  +  <f>o),  (121) 

V’2  =  B2{'ip-i  +  <^o))  (122) 

and  substituting  (122)  into  (121),  we  obtain 

{I  ~  B\  ■  B2)ipi  =  -Si(J  -|-  B2)(i>o-  (123) 

It  immediately  follows  from  the  combination  of  (117),  (112)  that  there  exists  a 
constant  Ci  >  0  such  that 

||G^l||oO  Cl  •  Halloo  ■  Pi  (124) 

for  all  p>  Q.  Combining  (124)  with  (116),  we  immediately  see  that  there  exist 
positive  constants  C2,  C3  such  that 

Pill  <  C2-P,  (125) 

11^1-^211  <  cz-p,  (126) 

for  all  p  >  D.  Therefore,  for  all  sufficiently  small  p, 

00 

{1-Bi-B2)-^=  ^iB,-B2r.  (127) 

m=0 

Combining  (127)  with  (123),  we  have 

=  (fliBi- B2)A  B^{I  +  B2)(I>o-  (128) 

\m=0  / 

Finally,  combining  (128)  with  (125),  (126),  we  obtain 

‘4^1  =  Bi{I  +  B2)(f>Q  +  0{p^),  (129) 

which  proves  (118).  The  substitution  of  (129)  into  (122)  yeilds  (119).  □ 
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Remark  2.35  For  any  u  €  the  function  v  €  C(R^)  defined  by  the  for¬ 

mula 

v  =  -{I  +  Gj)-^  ■Gru  =  Bj-u  (130) 

is  obviously  a  solution  of  the  Lippmann-Schwinger  equation 

(7  +  Gj)  ■v  =  -Gj-u  (131) 

on  the  domain  Q,j.  In  other  words,  for  each  of  j  =  1,2^  Bj  is  the  operator 
converting  the  incoming  field  into  the  scattered  field. 


Remark  2.36  Due  to  Remark  2.35,  (120)  can  be  interpreted  to  mean  that  the 
total  scattered  field  +  is  a,  superposition  of  five  scattered  fields 


(132) 

i=l,5 

where 

Vi  =  Bi-  <j)o, 

(133) 

V2  =  B2  •  4>o, 

(134) 

V3  =  B2  •  Bl  •  <f)Q, 

(135) 

V4  =  Bl-  B2-  <f>o, 

(136) 

V5  =  B2  •  Bl  •  B2  •  4>0’ 

(137) 

Each  of  the  five  scattered  fields  is  generated  in  a  scattering  process  described  as 

follows: 

1.  The  incoming  field  (f>o  gets  scattered  by  the  scatterer  in  fli,  generating  the 
scattered  field  v\.  We  denote  this  scattering  process  schematically  by  the  chart 

4>o  — >  — >  ui;  (138) 

2.  The  incoming  field  (f>o  gets  scattered  by  the  scatterer  in  Q2,  generating  a  scat¬ 
tered  field  V2 .  We  denote  this  process  by  the  chart 

(j)Q  - ^  0,2  - ^  '^2)  (139) 

3.  The  scattered  field  Vi  enters  the  scatterer  in  O2  o,s  an  incoming  field,  generating 
the  scattered  field  V4.  We  denote  this  process  by  the  chart 

4>o  — >  Oi  — Vi  — >  O2  — ^  V3’,  (140) 

4.  The  scattered  field  V2  enters  the  scatterer  in  Di  as  an  incoming  field,  and  is 
again  scattered  by  Oi,  generating  a  scattered  field  V3.  This  process  is  denoted  by 

(f)Q  — )•  rt2  — >■  V2  — >  rii  — »■  V4,  (141) 
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5.  The  scattered  field  V3  enters  the  scatterer  in  O2  os  an  incoming  field,  generating 
the  scattered  field  V5.  This  process  is  denoted  by  the  chart 

<f>o  — >  0,2  — ^  V2  — >  fli  — y  V4  — y  O2  — *  V5.  (142) 

The  following  lemma  is  a  restatement  of  Lemma  2.34  in  the  special  case  when 
0i  is  the  annulus  A{R,  h),  and  ^2  is  the  disk  D{R). 

Lemma  2.37  Suppose  that  R  >  0,  h  >  0  are  two  real  numbers  and  that  k  £  C'^ . 
Suppose  further  that  fo  :  D{R-\-h)  »-»•  C  is  an  incoming  field  to  the  disk  D{R-^h); 


in  other  words,  there  exists  a  £  Yk[R+h),  such  that 

■  Jkr  ■  oc  (143) 

inside  the  disk  D{R  +  h).  Finally,  suppose  that  ip  :  Rf  C  is  the  corresponding 
scattered  field  from  D(R  +  h).  Then  for  sufficiently  small  h  >  0  and  r  >  R-\-  h 

vfir,6)  +  0{h'^),  (144) 

i=i,5 

where 

vi{r,  9)  =  Ff''-  ■  Hkr  •  Spoc  -  a,  r  >  R+h,  (145) 

V2{r,  9)  =  Ff'^  ■  Hkr  ■  S^^k  •  r>R,  (146) 

^3(^5  ^)  —  Fg  •  Hkr  ■  Sf^  k  '  ‘S'oa  ■  <^5  r  >  R,  (147) 

^4(7’,  9'j  =  F g  ■  Hkr  ■  Spp  •  Sji  i^  •  cn,  r  ^  R  h,  (148) 

^5 (’'5  ^)  —  Fg  ■  Hkr  '  Sfi  f,  •  Sap  •  Sjip,  -  a,  r  >  R.  (149) 


The  following  lemma  is  analogous  to  Lemma  2.37.  it  is  a  restatement  of  Lemma 
2.34  in  the  special  case  when  Oi  is  the  annulus  A{R,  h),  and  ^2  is  the  disk  E{R). 

Lemma  2.38  Suppose  that  R,  h  are  two  positive  numbers  and  that  k  £  C"^. 
Suppose  further  that  fo  :  E{R)  1— >  C  is  an  incoming  field  to  the  scatterer  qE(R), 
and  that  ip  :  D(R)  i-+  C  is  the  corresponding  scattered  field.  Then 


i’iF0)  =  YMF0)  +  O(h^),  r<R,  (150) 

j=l 

with 

Ipi{r,  9)  =  Ff^  ■  Jkr  ■  Sap-  ffi  r  <  R,  (151) 

9)  =  Fg^  •  Jkr  •  r  <  R  +  h,  (152) 

^3(^5  d)  =  Fg  ^  •  Jkr  •  •  Spp  •  (3,  r  <  R  F  h,  (153) 

Ip  fir,  9^  =  Fg  ■  Jkr  ■  Sact  ■  ^R^-fi.k  '1^1  n  <.  R,  (154) 


9)  —  Fg  ■  Jkr  ■  ■  spa  ■  ■  fi,  r  <  R  F  h.  (155) 
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Remark  2.39  The  five  scattered  fields  xf^j,  j  =  1,5  in  (150)-(155)  are  generated 
in  the  scattering  process  described  as  follows: 


1)  <i>o^A{R,h)^rPi,  (156) 

2)  <^0  ^  EiR  +  h)-^  ^^2, 

3)  4>o  ^  A{R,  /i)  ^  E{R  Ah)^  V’s,  (158) 

4)  <i>o  E{R  +  /i)  — >■  'i/’2  A[R,  h)  (159) 


5)  (f>o-^E{R  +  h)^^2-^A{R,h)-^rp4^E{R+h)-^'il>5.  (160) 

3  Riccati  Equations  for  Scattering  Matrices 

It  turns  out  that  the  scattering  matrices  5"^.  and  5^^,  viewed  as  functions  of  r, 
satisfy  certain  matrix  Riccati  equations  with  respect  to  r.  In  this  section,  we 
derive  these  Riccati  equations. 


3.1  Scattering  Matrices  for  Thin  Annuli 

In  this  subsection,  we  obtain  approximate  expressions  for  the  four  submatrices 
Si30,  Sea  and  Sai3  of  the  matrix  (98)  when  the  annulus  A{R,h)  is  thin  (i.e., 
h  is  small).  In  this  case,  the  scattering  from  qA{R,h)  is  weak,  and  Spa,  Spp,  Sea 
and  Sap  assume  a  particularly  simple  form. 

Lemma  3.1  (Born  Approximation)  Suppose  that  R  >  0,  h  >  0  are  two  real 
numbers  and  that  k  ^  C'^ .  Suppose  further  that  <j>o  :  A{R,  h)  C  is  an  incoming 
field  to  the  annular  scatterer  qA(R,h),  and  that  So  is  the  restriction  of  So  on  the 
circle  Cr  C  RS'  ■  Finally,  suppose  that  S  •  R^  C  is  the  scattered  field.  Then 
for  small  h, 

S{r,e)  =  ^  ^  ^  —Ff^  •  Jkr  ■HkR-F-QR-So+  0{h^)  (161) 

for  all  r  <  Rj  and 

.H,.-F-QR-<t.§  +  0(h^)  (162) 

for  all  r  >  R-\-  h. 

Proof.  The  scattered  field  S  satisfies  the  Lippmann-Schwinger  equation  (see 
Remark  2.17) 

S{x)  +  k‘^f  Git(x,09(0t^(0<^^  = 

JA(R.h)  JA(R,h) 
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whoes  solution  can  be  approximated  by  the  sequence  {  0^?  fn  =  1,2,...  } 
generated  by  the  fixed-point  iteration  (72),  (73)  so  that 

*(l)  =  I  (164) 

JA(R,h) 

i’m+i  =  g-Gl-tpm,  m  =  l,2, ...;  (165) 

usually,  V’l  is  referred  to  as  the  Born  approximation  to  the  scattered  field  ip. 
According  to  Lemma  2.19,  there  exist  positive  numbers  e  >  0,  c  >  0  such  that 

\  ip(x)  -  ipi{x)  \<  c  ■  h'^  (166) 

for  all  h  <  e,  X  £  B?.  The  combination  of  (164),  (48)  and  (51)  yeilds 

j  fR-\-h  /•27r 

*(>-.«)  =  -.k'‘  /  ,(r',  «')*(/,«')>: 

4  jR  Jo 

oo 

Y,  (167) 

77l  =  — OO 

for  all  r  <  R,  and 

=  -k^  I  I  q{r',e')<f>o{r',9')  X 

4  JR  Jo 


oo 


Y  Jm{kr')Hm{kr)J^^^-^''>r'dr'de', 

m=~oo 

(168) 

for  all  r 

>  R  +  h.  Obviously,  for  any  function  g  £  ( 

U^[0,oo), 

r-R+h 

In 

g{t)dt  =  h  ■  {g{R)  +  0{h)), 

(169) 

and  combining  (167),  (168) 

with  (169),  we  obtain 

ipi{r,0) 

7  ^ 

=  -k^hR  Y 

(jm{kr)Hm{kR)  J\- 

q{R,6')(j>o{R,e')de'^^ 

je™' 

+Oih^). 

(170) 

for  all  r 

<  i?,  and 

ipi{r,d) 

=  '-k^hR  Y 

4 

[jm{kR)Hrn{kr)  fje- 

'^^'qiR,e')<f>o{R,o')de'^^ 

|e'”' 

(171) 

for  all  r 

>  Rl-  h.  Using  the 

i  notation  introduced  in 

Section  2.1,  (170),  (171)  can 

be  rewritten  in  the  form 

V’i(r,  9)  =  h^k^27rR  ■  Fg^  -Jkr-  HkR  ■  F 

■QR-<l>S  +  0{h^), 

(172) 
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for  all  r  <  R,  and 


Mr,0)  =  h'-k^2irR  ■  ■  hn  ■  H^r  ■  F  ■  Qr- 4>^  +  0{h\  (173) 

for  all  r  >  i?  +  ft.  Now,  the  lemma  follows  from  the  combination  of  (166),  (172) 
and  (173).  □ 

The  following  lemma  provides  the  desired  approximate  expressions  for  the  scat¬ 
tering  matrix  Sr  ^..  It  is  a  direct  result  of  formulae  (161)  and  (162). 

Lemma  3.2  Suppose  that  under  the  conditions  of  Lemma  3.1,  the  four  subma¬ 
trices  of  the  scattering  matrix  Spai  Spp,  Saa  Sap  defined  by  the  formulae 
(97),  (98)).  Then 


Spa 

1  TV 

=  hjeR-JkR-F-QR-F-^-JkR  +  0{h^), 

(174) 

Spp 

=  hjPR-JkR-F-QR-F-^-HkR  +  Oih^), 

(175) 

Saa 

=  hjeR-HkR-F-QR-F-^-JkR  +  0{h^), 

(176) 

Sap 

=  hj^R-HkR-F-QR-F-^-HkR-i-Oih^). 

(177) 

Proof.  We  will  prove  only  formula  (174),  the  proofs  for  the  rest  being  similar. 
According  to  Lemma  2.27,  an  incoming  field  fo  to  A{R,  ft)  has  the  form 


<i>0  —  Fg^  {JkrOl^  Hkr^^ ,  (178) 

with  a  G  Yk[R^h)i  fd  £  XkR.  Setting  /5  =  0,  and  substituting  (178)  into  (161),  we 
obtain  the  scattered  field  ^  outside  D{R  +  ft)  in  the  form 

^{r,9)  =  F^^HkrP  =  h-k'^R  ■  Ff^HkrJkRFQRF-^JkRcc  +  0(ft2),  (179) 

that  is, 

ITT 

IS  =  h-k^R  ■  JkRFQRF-^JkRa  +  (9(ft2).  (180) 

On  the  other  hand,  since  =  0,  (99)  assumes  the  form 

^  =  SpaOL,  (181) 

and  (174)  follows  from  the  combination  of  (180),  (181).  □ 
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3.2  A  Riccati  Equation  for  the  Scattering  Matrix  S 

In  this  subsection,  we  demonstrate  that  the  scattering  matrix  S~^  as  a  func¬ 
tion  of  r  is  a  solution  of  a  Riccati  equation.  To  this  end,  we  calculate  the 
difference5^^.;jji.  to  second  order  in  h,  by  observing  that  the  scattered 

field  from  the  scatterer  qn^R+h)  is  the  result  of  combined  effect  from  the  two 
scatterers  qD(R)  and  qA{R,h)- 

According  to  Lemma  2.20,  an  incoming  field  inside  the  disk  D{R-\-h)  has  the 
form 

<t>0  =  (182) 

with  Q  €  Yk[RAh)-  The  corresponding  scattered  field  outside  the  disk  D{R  +  h) 
has  the  form 

0  =  (183) 

with  €  Xk(RAh)-  By  the  definition  of  the  scattering  matrix  SRAh,k  (see  (87), 
(83)) 

^  =  ^R-k-h,k*^-  (184) 

The  following  result  is  an  immediate  consequence  of  Lemmas  2.37  and  3.2. 

Lemma  3.3  Suppose  that  R  >  0,  h  >  0  are  two  real  numbers,  that  (j>o  :  D{R  + 
h)  ^  C  is  an  incoming  field  to  the  disk  scatterer  in  D{R-\-h),  and  that  xj)  :  R? 

C  is  the  scattered  field.  Then 

^  =  {Sfik  +  h'-^P{JkR  +  S^^,HkR)FQRF-\HkRS^^,  +  J,r)}  a  +  Oih^). 

The  following  theorem  is  an  immediate  consequence  of  (184)  and  Lemma  3.3. 

Theorem  3.4  (Riccati  Equation  for  the  Scattering  Matrix  S~)  For  any  k  €  C'^ 
and  all  r  >  0,  the  scattering  matrix  Sfi^  :  Y^t  Xkr  is  a  solution  of  the  Riccati 
equation 

^  =  '-^k^iJkr  +  S;kHkr)FQrF-\HkrS;k  +  Jkr).  (185) 

dr  Z  ' 

3.3  A  Riccati  Equation  for  the  Scattering  Matrix  S'^ 

In  this  subsection,  we  derive  a  Riccati  equation  for  the  exterior  scattering  matrix 
Sf^k  defined  in  (90).  We  will  only  state  the  results,  since  their  proofs  are  quite 
similar  to  those  for  the  interior  scattering  matrix  in  the  preceding  subsection. 

According  to  Lemma  2.21,  to  E{R)  (the  exterior  of  the  disk  D{R)),  an  in¬ 
coming  field  assumes  the  form 

(f>o{r,  0)  =  Ff^Hkrfi,  (186) 
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with  ^  €  XkR.  The  corresponding  scattered  field  inside  the  disk  D{R)  assumes 
the  form 

rl>{r,e)  =  Fi^Jkra,  (187) 

with  a  €  YkR.  By  the  definition  of  (see  (90)), 

a  =  (188) 

Since 

E{R)=E{RFh)\JA{R,h),  (189) 

the  field  i/>  scattered  by  the  scatterer  qE{R)  can  be  viewed  as  the  field  scattered 
by  the  combination  of  scatterers  in  A{R,  h)  and  E{R-\-  h).  The  following  lemma 
is  analogous  to  Lemma  2.37,  and  is  an  immediate  consequence  of  Lemmas  2.38 
and  3.2. 

Lemma  3.5  Suppose  that  R,  h  are  two  positive  numbers  and  that  k  G  C'^ . 
Suppose  further  that  (j)o  :  E{R)  C  is  an  incoming  field  (186)  to  the  scatterer 
Qe{R))  that  ^  :  D{R)  i— >  C  is  the  corresponding  scattered  field  (187).  Then 

0  =  FQn^.,F-'  X 

{■Jk(R+h)SR^h^k  +  HkiR+H))}a  +  0{h^).  (190) 

The  following  theorem  is  an  immediate  consequence  of  (188)  and  Lemma  3.5. 

Theorem  3.6  (Riccati  Equation  for  the  Scattering  Matrix  S'^ )  For  any  k  G  C''*’ 
and  all  r  >  0,  the  exterior  scattering  matrix  Sf^.  :  Xkr  ^  Ykr  is  a  solution  of  the 
Riccati  equation 

-r^  =  +  SXMFQrF-\JkrSt,k  +  Hkr).  (191) 

dr  Z 
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